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ABSTRACT 


This report deals with the development of requirements for an approach 
control system and includes example applications to a jet transport air- 
craft. All of the techniques used have been known for some time, but the 
process of going from a list of guidance and control requirements to a 
system design has not been well documented in the past. 

The material presented is divided into two basic parts: a general 

discussion of approach control requirements, and a specific application 
resulting in the design of three alternative longitudinal controllers. 

The point of view taken is that the essential features of the system 
structure are the feedbacks themselves, their equalization, and their 
combinations to create control commands. Use is made of the fact that 
for successful systems the possible feedback structures are very limited. 
They derive primarily from guidance, control, and regulation demands; 
and secondarily from dynamic response characteristics desired by the 
pilot. From the systems view it is the satisfaction of these require- 
ments that is Important rather than the means employed. For this reason, 
most of the discussion in this report is equally applicable to automatic, 
manual, or hybrid manual/automatic approach systems. 
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SECTION I 


INTRODUCTION 

This report documents the development of requirements for an approach 
control system and includes example applications to a jet transport air- 
craft. Although all of the techniques used in this report have been known 
for sane time and in spite of the fact that there have been many pages of 
explanation devoted to design considerations, the process of going from a 
list of guidance and control requirements to a system design has not been 
well documented in the past. 

The material presented herein is divided into two basic parts: a 

general discussion of approach control requirements, and a specific appli- 
cation resulting in the design of three alternative longitudinal controllers. 
(These three controllers are compared and evaluated in terms of Category II 
approach success probabilities in Ref. 1.) 

A brief discussion is presented next to explain the point of view taken 
in this report. This will help orient the reader, and will also bring out 
some of the design ground rules . 

In approach and landing operations, the aircraft is but one element in 
a feedback control system. The essential features of the system structure 
are the feedbacks themselves, their equalization, and their combinations 
to create control commands . For successful systems, i.e., systems which 
demonstrate uniform, reliable, high quality approach and landing perfor- 
mance, the possible feedback structures are very limited. They derive 
primarily from guidance, control, and regulation demands; and secondarily 
from dynamic response characteristics desired by the pilot. From the sys- 
tems view it is the satisfaction of these requirements that is important 
rather than the means employed. In other words, the feedback loops closed 
is the central issue whether the closures are accomplished automatically 
or manually. For this reason, most of the discussion in this report is 
equally applicable to automatic, manual, or hybrid manual/automatic approach 
systems. Any differences come at a later stage when the feedback functions 
required are divided between animate and inanimate controllers and when the 
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subtle differences "between automatic controller and pilot dynamics are 
taken into account. 

Stated verbally, the key guidance and control requirements for low- 
level approach systems are: 

• To establish and maintain the aircraft on a 
specified spatial pathway or beam (e.g., 
localizer and glide path); 

• To reduce flight path errors to zero in a stable, 
well damped and rapidly responding manner; 

• To establish an equilibrium flight condition; 

• To limit the speed or angle of attack excursions 
from this established equilibrium flight 
condition. 

The regulation requirements are similar, i.e., 

• To maintain the established flight path in the 
presence of disturbances such as gusts, crosswinds, 
and wind shears; 

• To provide a degree of short-time attitude stability 
in the presence of disturbances. 

These requirements relate primarily to the relatively low frequency path 
modes of the air craft /control system. In essence, they define outer con- 
trol loops involving those vehicle motion quantities which define the 
desired equilibrium state of motion. More often than not, such outer 
loops, when closed about unmodified aircraft dynamics, do not result in 
stable, well-damped, rapidly responding systems. Instead, equalization 
of either a series or a parallel nature is needed to assist. Parallel 
equalization is most common and is achieved by the use of inner loops 
which feed back such quantities as attitude, angular velocity, and seme- 
times linear acceleration. These inner loops dominate the high frequency 
characteristics of the aircraft/controller system. 

To obtain a better appreciation of just what feedbacks the verbal 
requirement statements imply, we shall consider in Section II the deter- 
mination of feedback structures for a simplified lateral approach controller 
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and a more complex longitudinal controller. This -will be followed in 
Section III by a detailed development and analysis of three successively 
more complex longitudinal approach control systems. To make the discussion 
in Section III concrete, a DC -8-like aircraft will be assumed and numerical 
values will be used throughout. 
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SECTION II 


DEVELOPMENT QF FEEDBACK STRUCTURES 
A. LATERAL CONTROL 

A simplified lateral approach controller is shown in Fig. 1. The 
fundamental error signal in the system is the lateral displacement from 
the beam (y e ). This displacement from the beam is the difference between 
the beam*s lateral displacement (y c ) and the aircraft* s lateral displace- 
ment (y). The beam displacement (which is the commanded lateral displace- 
ment) is the sum of the desired lateral displacement (y- j_) and beam noise 
(n^). When lateral guidance is provided by a localizer, yj_ is the runway 
centerline, and thus equal to zero; for variable path systems y^ is a 
path command. Getting back to the error signal, y e can be converted to 
an angle (X) sensed by instruments in the aircraft via a relation involving 
the range from the localizer transmitter (i.e., without noise, y e = RX). 

It is noted that in Fig. 1 the measured lateral position error is 
contaminated by two kinds of noise. In addition to beam noise, there is 
receiver noise, n r . Typically, n^ is used to represent all unwanted inputs 
which are approximately stationary when represented as lengths (i.e., range- 
independent noise), and n r is used to represent noises which are approxi- 
mately stationary when represented as angles (range-dependent noise). 

Range -dependent noise includes the effective receiver noise, which tends 
to have a constant rms value at the output of the receiver and thus rep- 
resents a larger displacement at the longer ranges. An example of range- 
independent noise is main beam multipath transmissions. For the localizer 
these are caused primarily by fixed structures. For the glide slope, 
changes in ground reflection coefficients due to stratified wet and dry 
layers in the ground ("fixed" for a particular approach), and other devia- 
tions of the ground plane from an ideal reflecting surface are important 
causes. 

The receiving, filtering, gain changing, and other operations are 
represented in the transfer characteristic G\, which has the output 5^. 
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This is the effective command to the aircraft /flight control system. 

The latter comprises the airframe and associated inner-loop controllers. 

In the Fig. 1 diagram it is assumed that the aircraft controller 
feeds hack functions of hank angle, cp, and heading, t|/, to modify the 
basic airframe characteristics. The transfer characteristics &cp and 
G^ may he supplied hy the pilot and/or an automatic controller. To he 
explicit -we shall assume here that these functions are performed auto- 
matically. However, the properties subsequently developed for these 
transfer characteristics in order that guidance and control requirements 
he met are also incumbent upon the pilot if he is to play the same role. 

In the flight controller block diagram, note that a heading error, 
ijf e , is developed hy the insertion of a command or bias heading reference 
signal, \|f c . Note also that the relationship between heading and hank 
angle is given hy the simplified transfer function, g/U Q s. This 
simplification, as well as that between the flight path angle, 7 C , and 
heading, are consequences of assuming that the airplane is represented 
hy a three degree of freedom (spiral, roll subsidence) set of simplified 
equations of motion.* The total aircraft flight path angle, 7, is the 
sum of that commanded in the flight control system plus an increment due 
to crosswinds or gusts, Pg = Vg/U 0 . Finally this is converted into a 
lateral position hy multiplying hy U Q and integrating. 

For the present example we are concerned primarily with path modes; 
for these, the already simplified system of Fig. 1 can he further simpli- 
fied if we confine our attention to frequencies less than ocy Over this 
range of frequencies the flight path angle and heading are approximately 
equal; and, ordinarily, the amplitude ratio of the open-loop roll system 
(GqjO^a) is very large. With these simplifications the block diagram 
becomes that shown in Fig. 2. 


*These equations are the conventional three-degree-of-freedom lateral 
set with the simplifying assumption that (s — Yv)p is negligible relative 
to other side-acceleration quantities. See Section 6.7 of Ref. 2. 
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Figure 2. Simplified Lateral-Axis Block Diagram 


For this system the Equation of motion is 



V± + n b + 


G.r a 

R(t)n r +R(t) J t c + (- 


0 ) 


where p = d/dt. Because the range varies with time, this equation has 
time varying coefficients, although it is still linear. The range 
variation gives an increased sensitivity as the aircraft approaches 
the aiming point. Although a terminal controller which takes advantage 
of this time variation can he developed to give satisfactory accuracies 
at the end of the approach, the intrinsic time -varying gain must he offset 
hy a compensating time -varying gain jlf precision path control is to he 
maintained throughout the approach. This is done hy inserting a range- 
varying gain as one of the operations in such that ( 1 /R ) ( ) = G y , 
where Gy is a constant- coefficient operator. (Altitude or time, which 
are roughly equivalent to range when on a constant speed approach down 
a straight path, may he used in place of range.) Then the equation of 
motion becomes constant coefficient and can he written in Laplace 
transform notation as 

(~g s2 + ^ S + G y) y = G y (y i + % + Rn r) + + 0g) + Y % S Pg ^ 
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We shall study this equation in two ways. First, we will examine the 
steady-state characteristics in the presence of crosswinds to determine 
what is needed for windproofing. This shall be followed by consideration 
of the dynamics, using the characteristic equation, to determine the neces- 
sities imposed by path stability and response considerations. 

1 . Windproofing 

To represent a "nearly" steady crosswind, v c , in the steady state, 
assume Pg is a step function, v c /U 0 s. Then the steady-state characteris- 
tics of the system in response to this crosswind will be 



This can be made zero in several ways, i-or the second term to be zero the 
form sG^/Gy must have a net free s (or higher order) in the numerator. Common 
possibilities include 


sG<p, , sK cp 8 V ( V + 1) 

G y S " ^ S + K y ° r Ky s + K y 


(*0 


Addition of beam integration (i.e., adding a term K^/s to G^) provides another 
numerator s to offer further improvement. 

For the first term, if G^/Gy(0) is a constant, then a bias step command 
\|/ c /s must be introduced to the heading reference to just compensate for the 
crossvind. This requires either iterative trimming operations or precise 
knowledge of v c . A better technique is to provide a free s in the numera- 
tor of G^/Gy. Commonly used possibilities to achieve this include 
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K-s + K, 
i i 


K^y-S + Ky + Ky/s 


or 


% 


V 




s + K 


or 


G Jl . o 

G y 


(5) 


The type of windproofing actually selected depends to some extent on the 
feedbacks needed for dynamic control purposes, discussed next. 

2. Dynamic Requirements 

The effects of feedbacks on the dynamic characteristics of the system 
can be examined by considering the characteristic equation 

V 2 + V + s % = 0 ( 6 ) 

The low frequency forms for the controller transfer functions most con- 
ventionally used are 


G = 

9 T^s + 1 

^ = S ( 7 ) 

pv 

Gy = %s + % + -jL 


These forms are general and not all are used together. For example, the 
bank angle will not ordinarily be lagged if a proportional heading is 
present. The characteristic equation with these transfer functions 
inserted becomes 




(^ + S )] s2 + w) s + ^ <V V?> 


wSS . 


V 


( 8 ) 
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Consider, first, Eq. 8 with the path duration integral feedback, Ky, 
zero. The second-order system resulting represents an approximation to 
the dominant path mode of the system. For the beam to be followed at all, 
the constant term must always be present and have a positive value. For 
the path mode to have any damping the s term must also have a positive 
coefficient. As seen from Eq. 8 this can be provided by heading (K^), 
lagged bank angle (T^Ky), path rate (Ky), or by combinations thereof. In 
addition to providing damping, a heading feedback also provides attitude 
control that is inimical to mid-frequency windproofing. What happens, of 
course, is that an aircraft with a tight heading loop, when hit by a cross- 
wind, has a tendency to drift while maintaining a constant heading. This 
drift will ultimately be brought back by the Ky and Ky/s feedbacks, but 
only slowly. Path rate (y) on the other hand provides superior wind- 
proofing at the expense of heading. In the past y has been a difficult 
signal to obtain because of beam noise, so heading has been the typical 
path damping term. This situation is improving, however, due to the use 
of complementary filtering and the coming of scanning beams. 

Wow consider the complete third- order equation. The integral term, 
gKy/Kq^s, is present to assure steady-state windproofing if is not zero 
(see Eq. 5); it also suppresses steady-state lateral errors caused by 
crosswind shear. Ordinarily the first-order mode introduced by the inte- 
gral term has a very long time constant, given approximately by 


% 




Ky + TcpKy 


( 9 ) 


The corresponding approximate factors for the dominant second-order path 
modes are 
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In the above discussion we have considered only the path command or 
stiffening, path trimming, and path damping requirements, as these are 
fundamental to approach. In addition, there are requirements for attitude 
control and regulation. A bank angle feedback provides this function in 
roll. Also, near the touchdown point, tight heading control is needed to 
assure alignment of the aircraft wheel path with the runway to minimize 
landing gear sideloads. Just as the heading or path rate is required 
for damping of the path modes so is bank angle required as an inner loop 
for heading. Finally, to further Improve the total control and regula- 
tion precision, an extended flight control system bandwidth is desirable. 
This is achieved using roll rate feedback to the aileron. 


We assume in all of this, of course, that the yaw axis and, in 
particular, the dutch roll mode and any deleterious adverse yawing 
effects are taken care of by a suitable set of yaw damper and cross- 
feed loops. These will entail, in general, washed out yaw rate, side 
acceleration, and lag- lead aileron crossfeed (or their equivalent) fed 
to rudder. A block diagram for lateral control during approach incor- 
porating all of these features is shown in Fig. 3 . The path damping 
is provided by a combination of lagged bank angle and y derived from 
a so-called complementary filter. This appropriately mixes and filters 

a^- , cp, and, perhaps, a smoothed beam rate signal to obtain a broadband 

«y eg 

approximation to y. 


B. LONGITUDINAL CONTROL 


As another concrete example of system feedback selection we will 
consider a typical longitudinal approach control system.* The total 
system is shown in Fig. 4. There a distinction is drawn between a 
measuring subsystem and a control subsystem. The boundary is somewhat 


*It is noted in advance that the description of the longitudinal system 
will differ somewhat from that of the lateral system because the various 
simplifying assumptions used in the lateral case do not have longitudinal 
counterparts. Thus the longitudinal example will be more "involved" with 
algebraic detail, although the same kind of considerations (e.g., path 
damping, stability, windproofing, etc.) will still apply. 
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Measuring Subsystem 


Control Subsystem 



Figure lr . Measuring and Control Cystem for Vertical Plane Approach and Landing 




fuzzy and arbitrary in that certain other vehicle motion characteristics 
could just as well be included in the measuring subsystem portion. As 
shown, however, the measuring subsystem emphasizes the ground-to-air 
transmission of guidance data and the airborne decoding of these data into 
control system commands. 

The measuring subsystem has three important unwanted quantities. Two 
are the range -dependent and range- independent noises already described for 
the lateral example. The third unwanted input shown in Fig. b is the 
result of electromagnetic disturbances. These can come from a lead air- 
craft casting an electromagnetic shadow on following aircraft on the same 
approach beam pathway, from multipath transmissions of overflying craft 
(either direct or from side lobes), and so forth. 

As Fig- b emphasizes the measuring subsystem, so Fig- 5 provides 
a more detailed breakdown of the control system- Here the measuring sub- 
system is lumped into the beam smoothing and equalization block. Note that 
the hf signal, derived from a complementary filter combining barometric and 
inertial elements, comes from a different path than the deviation from the 
beam, d. While this h signal is shown in the controller in Fig- 4, its 
location in Fig- 5 indicates that it could just as well be in the measuring 
subsystem. 

In general, the development of an h-like signal can be accomplished 
using a combination of the beam rate signal with the outputs of baro and 
inertial elements in a complementary filtering scheme to reduce the effect 
of the beam noises on the derived rate. The inertial element can be as 
simple as an accelerometer or as complex as an inertial navigator. The 
latter are particularly appropriate for VTOL craft where the longitudinal 
position can be an important overall landing system loop- Inertial navi- 
gation equipment has also been proposed as a means to help reduce the effect 
of electromagnetic disturbances and noises (when properly combined with the 
other measuring system elements). 

As a simple example of complementary filtering, the composite signal 
derived from barometric rate of climb and accelerometer sensors, is shown 
in Fig. 6 . For simplicity the higher frequency lags inherent to the sys- 
tem are neglected. The actual rate of climb is h, and the various n's are 




Figure 5 . Block 














Figure 6. Derivation of Rate- of- Climb Signal 
unwanted signals and noises. The composite signal derived, h f , is given by 

n- + Tn 2 

= h + is + 1 ( 1 1 ^ 

Here, the major signal component is just h, while the noise terms are 
heavily filtered. The primary limitation on the time constant, T, is 
accelerometer noise, which is ordinarily very small at the frequencies 
of interest here. Consequently, the lag time constant, T, of the composite 
filter can be made quite large, with attendant reductions in noise. Thus 
the complementary filter can offer a potent means for the derivation of a 
good signal from a sum of signals which may be relatively poor individually. 
It is noted that the present glide slope beam is often too noisy near 
touchdown to permit its use for path rate computations, even with a com- 
plementary filter. This was recognized in the Fig. 5 block diagram, which 
has only baro and inertial elements called out. However, this situation 
may change as we gain more experience with Category II ILS, or with future 
systems. It is also noted that in addition to the noise situation there 
is a fundamental difference between the d and h signals. This is that d 
has a steady- state value of zero while h has a nonzero steady- state value. 
The significance of this difference will be presented later. 

Returning to describe the remainder of the control system block diagram 
(Fig. 5 ), it will be seen that both elevator and direct-lift control are 
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involved- For the sake of simplicity, a throttle control is not explicitly 
shown on Fig- 5 • As it stands, this block diagram is suitable for manual, 
automatic, or combined control because the mechanization of the several 
blocks is not specified. However, the control system is particularized to 
the point that attitude is fed back only to the elevator. The overall 
controller equations are given by 

S e Sp So 

5 e = % - (G d d + G 0 6 0) (12) 

Sf 8 -p 

S f = G d . di- G d d ( 13 ) 


Here, each G is a shorthand notation for the product of all the transfer 
functions in the blocks between the subscript variable and the superscript 
variable. For example, gS 9 is the product of the transfer functions for the 
blocks labeled: Beam Smoothing and Equalization, Elevator Input Equalization, 

and Elevator Equalization and Actuation — in short, everything between d d and 5 e 
Similarly, G d ? will be the product of the transfer functions for the 
blocks: Beam Smoothing and Equalization, Lift Control Equalization, and 

Lift Control Actuation. Using this notation the complete closed- loop 
approach system equations are given below. 


d = 


0 = 


+ C&g4 e ne + G^(N^ f + G^ e N§ e f f )n f 

[ G d| N & e + + Ge e N5eSf )]L + ^ 

sA + sGg e Ng e + G d e Ng e + G|f(4 f + G^ e | f ) 

W 

+ G§|( sNg e + G|f4J e )n e + Q§|( swg f + )n f 

^df ( sN 5 e + G d fN tfB e ) + G df( sN If + G d eN 8f5 e )] d i + + G d fN ^ f + G| e N®§ e )r} 


sA+sG 9 e N5 e + G d e »5 e + G |f ( f + G^N^ ) 


( 15 ) 


IT 


I I II I 


II II 


T) = Wg, U g 

A = Airframe- alone character is tie function 

0 d 

% e , e_tc - = Airframe-alone transfer function 

numerators 

• • 

N 5 e 8 fj etc. = Airframe-alone coupling numerators 

n e = Lumped noise effectively acting in 
elevator channel 

n^ = Lumped noise effectively acting in 
flap (DLC) channel 

Equations and 15 combine the controller equations with those of the 
vehicle, which is characterized by the transfer function numerators, 
coupling numerators, and characteristic function. These are summarized 
in Table 1 . Notice that the trim and atmospheric disturbances are 
denoted by a general disturbance input, r), and that the noises are 
lumped into an n e for elevator, and an np for the DLC channels respect- 
ively. Equations for other aircraft motion quantities, such as u, can 
be obtained from Eq. 15 simply by replacing the 0 superscripts in the 
numerator by the new variable. 

TABLE 1 

LONGITUDINAL EQUATIONS, TRANSFER FUNCTIONS, 

AND COUPLING NUMERATORS 


Equations of Motion 



d = -w + U 0 9 

h = -U 0 sin 0 O — w cos © Q + u sin © Q + U Q cos © o 0 
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TABIiE 1 ( C ont inued ) 


Characteristic Function (See note at end of Table) 



A A = 1 

Ba = -(Mq + Zw + W + Xu) 

C A = M q Z w - U 0 1V + X u (M q +Z w + U 0 M^) - + gM^ sin 0 O 

D a = -^(Z^-Uq^) - M u UoX w + M q Z u X w + g cos 0 o (I^ + Z u M^) 

+ g sin 0 o (M w -X u M^) 

E A = g cos © 0 (M W Z U -^Z W ) + g sin 0 o (M u X w — 


Numerators 


AgS 5 + BgS 2 



d 

B 


Ag = -z& 

Bg = z s(W + M q + X u) - X B Z u 

Cg = U 0 [M^(X 5 Z U -Z B X U ) + (ZgIV-ZwMg)] + M q (Z u X B -X u Z B ) + gMg sin 0 

B b = U 0 [X 5 (Z u M w -M u Z w ) - Z S (M W X U -X W M U ) - M^X^- Z^) ] 

+ g cos 0 o (Z u ]yi 6 -M u Z B ) + g sin 0 o (M u X B -X u M B ) 


® • «-y • * 0 • 

# = A^s" 5 + B d s 2 + C d s + D d 

"H ^ T] T] T) 

A d = +Z„ 

n ti 

B^ = -Z T] (U 0 M^+M q+ X u ) H-X^ 

= U 0 [M^(Z T) X U -X T1 Z U ) + (M^-Z^)] + M q (X u Z^ — ZyX^ ) - gM^ sin 0 O 

D d = U 0 [X q (M u Z w — Z-q^-) + Z T j(M w X u -X w M u ) + M^X-yZ-u — Z^jX^) ] 

+ g cos OoCZ^Mu-ZuM^) + g sin 0 Q (X u M 71 -MuX^) 
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TABLE 1 (Concluded) 

4 = A 5 s2 + B S s + C 5 


A® = Z 5 M^ + M 5 

B® = X 5 (M£Z u + M u ) + Zg^-M^) - Mg^+Xj 

C S = ^(^u-Mu^) + Z&^-MwXu) + M B (Z w X U -Xw Z u) 

K = A ®S 2 + B ® s + C ? 

A ® = -(Zr^ + M,,) 

B® = + [-X T1 (I^Z U +M U ) - Z^M^-M^Xu) + M^ + X*)] 

c ® = -Vfr^-VW-W-VzA-xA) 


Coupling Numerators 


„0 d 0 d 

N S e S f " A 5 e 8 f s 


+ B, 


0 d 
'8 e 5f 


A 8 e Sf = ( z 8e M Bf -Z Sf M &e^ 

B 5 e 6f = ^(Zg^-Zs^) - Z^^-Xs^) +M u (X 6f Z 8e -X 8e Z 5f ) 

<1 = A*g®s + B*® 

A* f = Z^Mg - M^Zb 

B* g = -XjZ^Mg-M^Zg) + Z^Mg-M^) - M u ( X^Zg - Z^Xg ) 


Note that some of the transfer function numerators are defined in terms 
of error rate (rather than error) in order to avoid the confusion that can 
arise from transfer function numerators having denominators. Thus, the 

transfer function for d/S 0 is defined to he N^/A^^g, and the transfer func- 
tion for d/5 e is N§ e /sAi on g* By not using the symbol N^, all numerators 
(Ng^) remain "pure” polynomials (rather than ratios of polynomials). 
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Having described the control system block diagram (Fig. 5), and 
defined the airframe transfer functions (Table 1 ), it is now appropriate 
to give a "verbal" description of the fundamental guidance and control 
requirements for the longitudinal approach system. 

• Guidance Requirements 

The aircraft must follow the beam commands 

• Regulation Requirements 

• The aircraft should be maintained close to 
the beam in the presence of winds, gust dis- 
turbances, internal biases in the equipment, 
measuring system noise, etc. 

• Aircraft attitude should be kept stable and 
"solid", i.e., relatively constant, in the 
presence of the disturbance environment. 

• Implied (Mechanization-based) Requirement 

Elevator must be used for trim adjustments (i.e., 

DLC is not to be saturated due to trim changes) 


Most of these qualitative requirements can readily be translated into 
required feedbacks by considering certain steady-state aspects of the 
system equations, followed by simple stability and response arguments. 
We shall consider the steady-state features first. 


Guidance Requirement 


The most fundamental guidance requirement is that the 
aircraft acquire the beam when the system is engaged. In 
other words, the deviation from the beam, d, must ultimately 
become zero when the system "input" is an initial condition, 
d( 0 +). The response transform of an nth order system with 
characteristic function 


sys 


s n + 


a-j s 


n-1 


a n-1 


s + a. 


n 


(16) 


to an initial condition of position, d(0+), is readily shown 
to be 


(s 31 " 1 + a-jS 11 " 2 + ... + a n _ 1 )d(0+) 
s n + a-|S n-1 + ... + a^s + a n 


( 17 ) 
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Assuming that the final value theorem holds , 


lam 
t — 


a(t) 


lim 

s-*^0 


sd(s) 


11m 
s — *~0 


(s n ~^ + a^s n "^ + ••• + a n _ <1 )d(0+) 
s n + a-j s n ^ + ... + a n _'| s + a n 


08 ) 


which will be zero only when / 0. Consequently the charac- 
teristic function must contain a constant term for the system to 
acquire the beam. The characteristic function is 


A sys ~ ^long + 


sG^| e + G^4 e + G d f 04 f + G^ e f f ) 


( 19 ) 


0 d 0 d 

As s approaches zero, sA long , sNg e , Ng, and Ng e g f approach sE A , 
sCg £ , D^, and Bg e g f , respectively. So, as s approaches zero in 
Eq. 19, 


im A 

s -*-0 sys 


3^0 { s [% + G 6 e ( s ) c 6e] + G d e ( s ) D 6, 


+ G* f (s) 


d 


D 6f + G^B5 e5f 


( 20 ) 


For the attitude control to be significant in this expression 
attitude feedback, Gg e , would have to contain an integral term. 

This would conflict with the desire to drive d to zero for other 
inputs. Consequently G| e (0) will be made either a constant or zero. 
Also, for reasons which will be described later, the DLC control 
of path deviation, G^*, should have one more free s than the 
elevator/path deviation transfer function, G^ e . This leaves 
the G^ e (s)D^ term, which will satisfy our need for a constant in 
A S y S if Gge Contains a proportional (K^ e ) term. Thus, all of this 
justifies the intuitively obvious requirement for a proportional 
feedback of path deviation to achieve static stability relative 
to the beam. 

Another guidance requirement is that the system follow guidance 
commands, dj_. These might arise from a one- step beam scheme or 
even a higher order curvature command in more advanced systems. 
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In following guidance commands, d^, the error d e is given 


by 


d £ = d ± - d 


sA + sG| e N& e + ]^ e (4 e ~ ' jji ) + ( W bf + G 6 eN 5eSf ) (ffi ~ G bj) 

SA + sG^ e N| e + G| e 4 e + <fi f (uff + ^bjs t ) 

o c; 

Letting G^ = G^ + for "both S e and then gives., 

sA + sGg e N5 e + Y|^ e + + G| e N| e f f ) ^ 

sA + sGg e Ng e + G| e 4 e + G| f (4 f + G^ e N|J f ) 


d e = 


d ± (21) 


( 22 ) 


If the commanded path is given hy a power series in time, i.e.. 


d i (t) = d 1 + dgt + d^t 2 + • 


then 


di do 

d i(=) - -r + -i + -r 

s 


2d^ 


+ . ■ 


(23) 

(24) 


and the lowest order term in s will “be (n— 1 ) Id^/s 11 . Using the 
final value theorem, 


lim 

^-steady state ~ s-*-0 


(sA+sG^^ + Y^N^ + llfl 

( N 5f + G 9 eN &eBf)l 

'(n— 1 )’dn 

jsA + sa| e w| e + G^ e 4 e + G| f < 


s*' 1 

1 


(25) 


For the system to follow a command path of nth order the 
numerator of the curly bracketed expression in Eq. 25 must 
contain a free s 11 . From the sketch it is apparent that a sys- 
tem which is stabilized on the first segment of a two segment 
glide path system must follow a ramp function in d without 
steady-state error if it is to successfully transition from 
path 1 to 2. So a free s 2 is needed in the curly bracketed 
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portion of Eq. 25. Using the same limiting properties as in 
the beam acquisition case. 


s 



11m 

s— 



+ G° e 


(s)C? e ] 


+ X 


5 (s)d£ 


&e( 

d ^~ /_ o e 


Gd e (s)D5 P 


( 26 ) 


5 

typically is either zero or contains a single free s, so 
the total numerator in Eq. 2 6 has a net free s. Then, to 
provide the second numerator free s (needed to satisfy the 
steady- state requirement) the path deviation/elevator trans- 
fer function must have an integral term, Krj e /s. This is 
also obvious intuitively from examination of Fig. 5 . There 
it is plain that a steady- state signal must be developed at 
the d c point to offset a change in the steady- state output 
of the Baro- inertial Smoothing and Equalization block. 


Regulation Requirements 


The path deviation response to an external disturbance 7 ] 
is given by 


N d + Gg e N d £ 

Ti y T] Oe 


sys 


(27) 


where T| may be a w g or Ug wind disturbance; and either can have a 
constant component, i.e., t](s) = T]^/s + ... . For the constant 
component to have no long term effect requires 


2k 



( 28 ) 


'ss 


. (n~ + Gg e N^ f 

s— ►O 3 a < T)(s) 

( ^sys 


= T] 


lim 
1 s-*-0 


= 0 


This will he zero when there is at least one net free s in 
the numerator. 


Case 1 : tj 1 = w g 


- L \r t 


Sss 


= w, 


11m 


S s— 0 


(29) 


r^e^d 

G d D 5 e 


Because the numerator is a constant as s— *-0, the 
free s needed to make the whole thing approach zero 
must come from the denominator. Consequently, a Krj e /s 
component is needed in the G^ e control path. A similar 
argument applies for trim changes resulting in Z c and Mq 
lift and pitching accelerations applied to the vehicle. 
Also, note that any shear component to w g will result in 
a steady- state error even with an integral controller. 

Case 2 : = u g 



u g 


lim 

s-*-0 


s 


I c ug + G 6 eA u g 5e 


(50) 


A free s occurs naturally in this numerator, so for d u 


to he zero requires only that A svs (0) he a constant or B 

- 3 _ mu • — /-tries f ~ \ 1 ^ 7y~c\ ■ 


?SS 


lower order in s. This means that G; 
satisfactory. However, head or tail 

t 


winds invar iabl; 


is 

.y have 


a shear component, so a K^ e /s is desirable to help counter 
these. 
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These elementary considerations indicate that the most undesirable 
disturbance inputs are wind shears. In principle, the worst of these is 
a shear normal to the flight path, for this will cause a steady- state 
error in path. In practice, however, all shears occurring near the 
terminal condition give cause for concern because their effects are 
countered primarily by the integral control. This is inherently slow 
in action, as will be appreciated better with the aid of the concrete 
examples of the next section. The promise of Eero steady- state error 
is, accordingly, more academic than real. 

Additional "steady- state" requirements on regulation against disturbances 
can be derived if a shorter time scale is presumed. This can be done by 
considering the two-degree- of -freedom short-period characteristics instead 
of the complete three- degree- of- freedom equations. Any "final" values 
found using the short-period equations apply for time intervals which are 
large compared with the system* s settling time, but not so long as to be 
comparable with phugoid periods. Thus the short-period approximation is 
valuable to treat some mid- frequency response properties. 

Using the short-period approximation (see Table 2), 



This equation is most pertinent for a w-gust disturbance,, When = w g , 


^Wn 


=ss 


= w g s- 


1m | V [ s ~ ^ + %)] + G 8 e ^e^- ^ 


■^sys 


(32) 


A S y S (0) will be a constant if K^ e is the lowest frequency feedback, 
or if the time span considered is relatively short, such that a K^ e /s 
control will have little effect. Then, to obtain a net free s in the 
numerator, Gg e must, itself, contain one. Thus we can establish a desire 
for pitching velocity rather than attitude feedback to improve the gust 
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TABLE 2 

SHORT PERIOD EQUATIONS AND TRANSFER FUNCTIONS 


Equations of Motion 

[(l-Z^s-Zjsd + 

2 

(M^s+M w )sd + [s — 


Z 8 
a 

(M q +M i )s-%]9 



0c)»e-Z6f»f- V 


«6e 5 e 


“ M Ti 

■n 


A B e — 6 f direct crossfeed and a possible a z feedback to give 

the Z^* are taken into account in these vehicle equations; and = 0 . 


Characteristic Function 


A 

Numerator* 


4 




2 ) 2 


= -Z. s — (M +M- 


N d 

*n 


K 


.0 

5 “ 
n g = 


Si)' - (*t-^ Z a) 

z J s2 - < M i +M ;) s -( , u-^ z < I ) 

s |[l%(l-Z£) + Z 5 M^]s - (M & Z w -Z 5 M w ) 

-s jtM^I-Zh) + Z^s - ( Vr - Vr )) 


Coupling Numerators 


N. 


_9 d 
S e Sf 


N^ - = 

T| 5 e 


-ZsfMBe 

*%e Z >l - 


Modal Response Ratio 


PV1 


L Zw 

\ + G»f 5* 

1 1 W 


V 

J d 


si 


*For these short-period equations it is possible to define N§ without 
having to resort to numerators that are ratios of polynomials . 
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regulation at mid- frequencies . This permits the drift from the "beam 
caused by a normal gust to be reduced by virtue of the aircraft* s -weather- 
cocking tendency. This would not be possible if the aircraft rigidly 
maintained its pitch attitude. 


Implied Requirement on 5f 


The direct lift control, be it spoiler or flap, can have 
only a very limited control power compared with the elevator- 
wing combination. Consequently the longer time (near steady- 
state) control should be elevator to avoid saturating the 
DLC . The time scale of interest here is relatively short, 
so the short-period equations (Table 2) can again be used 
to define a kind of "short time steady-state." For a trim 
change defined by an incremental lift, Z^, and pitching. 

Mo, accelerations, the flap deflection will be, 


6 f .= _G d‘ 
= 


Sf d 


- [s 2 — (M* + lfe) s — Ik] - Ge e M 5 e j z 0 + (~z a + G^ e )M Q 


A, 


sys 


(33) 


.Se 


If G 0 C is presumed to have a free s, and Z D and M Q have 
constant terms, e.g., Z Q (s) = Z 1 /s + ..., then (using Eq. 33 
and the short -period version of Eq. 19 ) 


6 f 


ss 


la |4 f (Ma Z 1 - Z a Mi) 

s — 0 






(3k) 


For 5f ss to be zero then requires Gj f to be one order higher 
in s than (§e. 
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From the above steady- state (long time as well as short time) con- 


siderations we have shown that the minimum forms of the feedbacks needed 
to satisfy reasonable steady- state guidance, control, and regulation 
requirements are: 


li 

<D 

Be 

S a. 

(55) 


G e e * 

V 

(56) 

and, if DLC is used, 

4 f - 

K^jf or higher terms in s 

(57) 


Having disposed of the feedback requirements imposed by steady-state 
considerations, we shall now turn to a short discussion of the higher 
frequency feedbacks. The simplest of these is the attitude transfer 
function, Gq . When the short term attitude regulation requirement is 
considered — implying attitude stiffening at short-period frequencies 
to provide a craft that is stable and solid in attitude — the Gg e of 
Eq. 36 is modified to a simple washout. This retains the free s in the 
numerator of Gg e , with its favorable consequences in regulating against 
wind disturbances, while still providing an attitude feedback at short 
period frequencies. The washout time constant 1 /T w must, of course, be 
such that l/T w is less than co S p. Further, taking higher frequency effects 
into account, a pitch rate feedback Kgs would be desirable to provide 
greater short-period attitude damping and attitude loop bandwidth. This 
would then permit improved altitude loop gain margins and, thereby, a 
greater altitude bandwidth. Consequently the desired general form for 
the pitch attitude feedback is 



v 

a + (1/T W ) 


+ K-s 


(58) 
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The remaining requirements for the general form of the feedback 
controller dynamics can be derived by analogy with the lateral case. 
Thus, with the attitude feedback washed out at low frequencies, 
feedbacks are needed to Improve the path damping. As is evident from 
Table 3, which shows the closed-loop characteristic function coeffi- 
cients, this can be a useful feedback to both elevator and DLC. In 
fact, because the high frequency limiting factors on DLC and elevator 
closures can be somewhat different, a system using the d feedbacks to 
both controls is desirable. Recall, however, the implied requirement 
that G^ f be one order of s higher than G| e . Thus, appropriate general 
forms for the path deviation feedbacks are 



5 e 

% 


K a e 


K; 


(39) 





s 


(^ 0 ) 


Without going further into all the ramifications and justifications 
behind these selections, a general summary of the feedbacks, their pur- 
poses, and qualitative requirements is provided in Table 4. 
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TABUS 5 


CLOSED-LOOP CHARACTERISTIC FUNCTION COEFFICIENTS 
USING SHORT-PERIOD AIRCRAFT EQUATIONS 


. 5 , 


0 


4 

s : 


s 5 ; 


% o -z£> - <i -Sk>(v>fc+T%) + K 0 [ V’ - z h> - + vl f ) 


d-^) 


L \ 

T* ’-d 




y , r\/r i, - ! hJ K ® , - z i' 


- Z8, 


4' + (^-« a -«i)K|< 


- Z. 


Of 


K 4 f+ (i-« 4 - M a)*r 


- K i 




Z M 


(1 - z h‘){{- t (VT^) " ( K e + S) )4 e + K d e (f ' 

1 / ^6f* 

! /-h/r i n/ . N i - 1 - 


■y \ e V 


1 


gr (V*fe) + 


ZSfMSg 


1 - ZI 


2h 


Sf / K 0 ' 


tfe-\ 


- M-] - 
a/ 


fe 

d 


J- (M +M-) + M - vr 2 Z 

T w q- a \ a z 5 e a 


M r 


Mb ^ 


+ Z 5f 4 f 


M6f 
5 f 




+ 4 f i [*.-% la) 


o 

s : 


z M K f 


(M +M 


T w v q. 


[.) + (m - 

a ' la 


Z 5 e “ 


+ K*?e 1 


IM - 

\ a 




Of i L ®^f 


<*’ M“ a z b 6 “ a /j + z& f r d z 5f z 7 


r ’= 


Mb, 


■(M„ +M-) + 
. q a' 




TABLE 4 


LONGITUDINAL FEEDBACKS — PURPOSES AND QUALITATIVE REQUIREMENTS 


FEEDBACK TERMS 

FUNCTION 

REQUIREMENT 

REMARKS 

Ke 

Short-Period Attitude 
Stiffness 

9 in short-period frequency 

range (l/T w < co sp ) 

Windproofing and attitude 
stiffening conflict 

K e 

Short- Period Damping; Path 
Loop Bandwidth Extension 
Capability 

0 in short-period frequency 

range 

T 

■ L w 

Short Term w Windproofing 

e> 

s[s — (M : +%)] + G^M 5 e 
small over w g frequencies 

K a 

Path Acquisition and 
Stiffness 

Windproofing (u g step, 
Wg pulse) 

sets dominant path mode fre- 
quency; made as large as 
possible consistent with 
stability and limiting 

Altitude control bandwidth with 
d,e -»B e is limited by 1 /Tgp; 
can be increased significantly 
by d,d, or h Therefore 

5f can be big help in approach 
and flare precision. 

I'd 

k 

J 

Higher-Order Path Following : 
Trim 

Windproofing (w step) 

o 

K d/ K d se ' ts trim response time 

5- 

Kg 1 = 0 for B e to handle trim 

K d or K h 

Path Damping 

Sets path mode damping ratio 

0.64 < £ . . - < 1 

13 path mode 

d suffers from beam noise; h 
requires trim bias to offset 
steady- state sink rate. 

K e 

0 in long- period frequency 

range 

Conflicts with windproofing. 




SECTION III 


DEVELOPMENT AND ANALYSIS OF EXAMPLE SYSTEM3 


To provide concrete examples of approach systems, three longitudinal 
approach controllers shall be developed in this section. (The results are 
used in Ref. 1 to compute average performance measures and probabilities of 
approach success . ) The three systems developed can be considered as competi- 
tors throughout the analysis; that is, they can be considered in terms of: 
dynamic characteristics, performance measures and probability of approach 
success . 

The aircraft to be controlled will be a DC-8 defined by the landing 
approach configuration parameters given in Table 5 . The aircraft transfer 

TABLE 5 


DC -8 PARAMETERS FOR LANDING APPROACH CONFIGURATION 



geometry and 


LONGITUDINAL 


LATERAL 

1 INERTIAL PROPERTIES 


STABILITY AXES 


BODY AXES 

h 

(ft) 

0 

*u 

( 1 /sec) 

-0.0575 

Y v 

(l/sec) 

- 0 .0887 

M 

(-) 

.204 


( 1 /sec) 

0.156 


( 1 /sec) 

0 

VTo 

(ft/sec) 

228. 

X *e 

(ft/sec 2 /rad) 

0 

Y &r 

( 1 /sec) 

0.051 

?o 

(deg) 

-2.8° 


(ft/sec 2 /^) 

0.106 


( 1 /sec 2 ) 

- 1 .4o 


(lb/ft 2 ) 

61.8 

Zu 

( l/sec) 

-0.285 

h 

( 1 /sec) 

- 1 .04 

S 

(ft 2 ) 

2758 . 

Z* 

( l/sec) 

^•750 

L r 

( l/sec) 

0.474 

b 

(ft) 

Ilf-2, 4 

Z w 

(-) 

0 

L &a 

( l/sec 2 ) 

1.15 

c 

(ft) 

22.16 

z 5e 

(ft/sec 2 /rad) 

- 9.25 

L 6 r 

(l/sec 2 ) 

0.159 

W 

(lb) 

180,000. 

Zgrjl 

(ft/sec 2 /$) 

- 0.00097 


( l/sec 2 ) 

0.568 

m 

( slugs ) 

5,580. 

M a 

( 1 /sec- ft) 

0 

N P 

(l/sec) 

-O .029 


(slug- ft 2 ) 

5.2 XlO 6 


( 1 /sec-ft) 

-0.00461 

N r 

(l/sec) 

-0.257 

*7 

( slug— ft 2 ) 

5.8 XlO 6 
£ 


(1/ft) 

- 0 .00085 


( l/sec 2 ) 

0 

I Z 

( slug— ft 2 ) 

6.6 x 10 


(l/sec) 

-0.594 

w 8r 

(l/sec 2 ) 

-0.568 

Y xz 

(slug-ft 2 ) 

0 

M8 e 

( l/sec 2 ) 

- 0.925 




*CG 

(# c) 

25.2 


(1/sec 2 /$) 

0.000625 




&F C 

(deg) 

50 

% 

(1/sec 2 ) 

- 1 .05 




a o 

(deg) 

0.62 


(l/sec) 

-0.1956 
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function characteristics for control inputs are shown in the Bode plots 
of Fig. 7. The notation in the numerical transfer functions shown on 
these plots is a shorthand in which (l/T) represents a first-order term 
with time constant T, and [£, o^] represents a second-order factor with 
damping ratio £ and undamped natural frequency o^. For example., in 
Fig. 7 , 1 /T 01 = 0.101, i/Tq 2 = 0.646, £ sp = 0.626 and ao sp = 1.251. The 
same notation is used in the complete compilation of transfer functions 
and coupling numerators given in Table 6 . 

The control equations and functions accomplished by the three systems 
are given in Table 7 . The systems are arranged from "A" to "C" in order 
of decreasing complexity and capability with "A" also standing for 
"advanced” and "C" for "conventional." All of the systems can acquire 
and maintain position on a straight line glide slope beam with well-damped 
path mode responses. However, System C is not suitable for following 
higher order paths with zero steady state error. The major distinction 
between Systems B and C is in the 6 feedback, which is washed out at very 
low frequencies on System B and not at all on System C. The washout is 
intended to improve the w_ windproofing, the steady-state following of 
higher order paths, and to remove the effects of any steady- state 0 
biases. This is achieved at the expense of a slight amount of path 
damping and bandwidth. Consequently, the superiority of System B in 
Wg windproofing and steady- state operations may be offset, for other 
inputs, by its smaller bandwidth. 

System A is representative of the elevator axis of an advanced 
controller, typical of the forthcoming generation of low level approach 
and automatic landing systems. Appropriate feedbacks exist for all the 
functions listed for longitudinal control in Table 7 . 

All of the example systems can be improved by the addition of an airspeed 
control. However, as will become evident later, the airspeed control proper- 
ties of the three systems as they stand are nearly identical. Further, this 
similarity will not be changed if the same airspeed controller is added to 
all systems. Consequently, for the sake of simplicity, we have not provided 
airspeed control loops. 

The control equations for the three systems are given in the last row 
of Table 7 in terms of This is the commanded elevator which must 


3 ^ 







TABLE 6 


LONGITUDINAL TRANSFER FUNCTIONS FOR THE DC-8 IN LANDING APEROACH CONFIGURATION 

Denominator 

A = [O.IO , 0.167][0.626 , 1 .231] 

Numerators 

S e Control Input 

Ng e = -1 .258(4. 03) (-4. 082) 

N5 e = -9. 25(23. 34)[o. 107 , 0.198] 

Ng e = -0.9151(0.101 )(0. 646) 

N& e = 9-239(0.042)(-3.607)(4.397) 

Ng e = 9. 25(0. 035) (-3. 606) (4. 396) 

u Gust Input w Gust Input 

—2 _£ 

< = 0.0373(1. 543) [0.599, 0.857] < = -0 .i 360(0)[0.486 , 0.795] 

§ & 

Nu_ = 0.283(0)(0)(0.594) hJ = 0. 3498(4. 379)[0.1 18 , O.166] 

n£ = -0.0002406(0) (5. 424) N^ = -0.001755(0)(-1 .475)(-0. 000464) 

n£ = -0. 2845(0. 007)[0. 386 , 1 .027] n| = -0.7425(1 .243)[0.08l , 0.2i4] 

U g g 

Nu = -0.283(0)[0.384 , 1 .025] nJ = - 0 . 75(i.234)[0. 091 , 0.2i4] 

g & 

Coupling Numerators 

N? e u g = -0.03413(1.751) 

N 5 e u g = -0*2612(0) 

N 5 e u g = 0.345(1 .102)(22. 243) 

<£ g = 0.1223(0) 

NgeWg = -O. 01624(0. 092)(39*955) 

»5eWg = -i.258(-0.696)(23.626) 
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TABLE 7 


CONTROL EQUATIONS AND FUNCTIONS PERFORMED BY THE EXAMPLE AUTOMATIC LONGITUDINAL SYSTEMS 


FUNCTION 

SYSTEM 

A 

B 

C 

Attitude Control 
and Regulation 

Short-Period Attitude 
Stiffness 

0 B e i 31 short-period frequency range 

Short- Period Damping; Path 
Loop Bandwidth Extension 
Capability 

0 -*-S e in short-period 
frequency range 



Short Term Wg Windproofing 

High-frequency 0 washout 

Low-frequency 0 washout 

No 0 washout 

Path Control and 
Regulation 

Higher- order Path 
Following 

Trim 

Windproofing (wg step) 

Jd dt S e 

Low-frequency 0 washout 


Path Acquisition and 
Stiffhess 

Windproofing (ug step, 
Wg pulse) 

d — & e 

Path Damping 

h 5 e 

0 5 e In long-period 

frequency range 

e — 6 e 

Control Equations 

+ K d s 

~Sp = — -d p + Kch 

C s(T f s + l) e h 

K e s 

+ e + k. e 

3+i /i„ 0 

K,d K ft s 

t - d e , 6 _ 

-Se c T f s + 1 s + 1/T^ 

K d d 

_5 e c ~ T f s + 1 + K 0 e 




be converted to actual elevator deflection by the flight control and 
surface actuation systems. The actuator characteristics (Y a ) will be 
approximated here by a first-order lag with 1 /T a equal to 15 sec 1 . 

Although grossly overs amplified, this is an adequate approximation to 
the actuator properties for the low frequency range of primary interest. 

Block diagrams corresponding to the control equations are given in 
Figs. 8 and 9 for Systems A and C. A block diagram for System B would 
be essentially the same as that in Fig. 9 , with the replacement of the 
attitude feedback, Kq, by the transfer function, K 0 ^wo s / ( ^wo s + 1 ) • 

Having defined the systems qualitatively, we shall now turn to the 
quantitative descriptions which will be developed for each system. 

First, the closed-loop dynamic characteristics will be considered In 
terms of Jco Bode and Bode root locus plots for the systems as adjusted. 

Then, typical time histories of the systems will be discussed for command 
and disturbance inputs. 

Conventional System (System C) 

The analysis of System C will be described in two steps: the closing 

of an attitude inner loop followed by the closing of a path deviation outer 
loop. The connections and interplay between these two operations is central 
to the synthesis procedure, because the attitude loop provides the equaliza- 
tion necessary for the path deviation loop to be closed such that rapid, 
stable, well-damped responses result. Consequently, we will take some 
pains to point out these connections in a discussion of the analysis. 

A jco-Bode plot for the open-loop attitude to elevator transfer function 
and a closed- loop Bode root-locus for the same system are shown in Fig. 10. 
The Bode root-locus comprises both real and complex roots, the real roots 
being shown with the heavy line (so-called siggy Bode plots), whereas the 
complex roots are dotted. (The values of real roots and closed- loop undamped 
natural frequencies are read using the abscissa, while the closed- loop 
damping ratios are parameters along the complex branches. Gain, of course, 
is the ordinate.) The Bode root-locus shows that as the gain, K 0 , is 
increased, the damping ratio and total damping of the phugoid mode are 
also increased, resulting (at high gain) in two real roots which approach 
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Figure 8. Block Diagram of Advanced System [A] 










TR- 182-3 







w,a, Isl log scale 


(rad/sec) 


Figure 10. System Survey for Attitude Loop Closure for System C 





the zeros, — i/Tq^ and — 1 /T g^. For the same gain variation, the closed-loop 
short-period damping and damping ratio are decreased. Thus, the increased 
phugoid damping is obtained at the expense of the short-period damping. 

For relatively low gain, the decrease in short-period damping is reflected 
primarily in the short-period damping ratio because the short-period 
undamped natural frequency is essentially unchanged. At moderate and 
high gains, the short-period natural frequency increases, resulting in 
an even faster decrease of short-period damping ratio with gain. Finally, 
the closed- loop actuator characteristics are but slightly modified from 
the open-loop properties for reasonable values of gain. For the nominal 
zero dB line shown (K Q is —3.65) the phugoid roots are shown as 1 /T^ 
and 1 / T 02 (^ e prime indicating one loop has been closed and the Tg^ 
and Tg^ indicating the zeros being approached), and the closed- loop 
short-period damping ratio is 0.1 84 . With this closure, there is a 
wide spread between i/Tq^ and 1 /Tq 2 , yet the short-period damping ratio 
is still moderate (£ S p = 0.1 84 gives 0.6l cycles to half amplitude). 

Only the poles of the glide slope deviation to elevator transfer 
function are modified by the closure of the 0 to be loop; the numerator 
is still that of the airframe alone. The deviation to deviation- error 
open- loop transfer function, with the attitude loop closed, is shown in 
the system survey in Fig. 11, where the jco-Bode reflects the open-loop 
zeros of the numerator, Wg e (s), the closed- loop poles resulting from 
the attitude closure, and a beam noise smoothing filter (Yf) with a 
time constant of 0.5 sec. The Bode root-locus of this (outer) loop 
is also provided in the figure. Here it is seen that the very- low- 
frequency mode stemming from the free s advances toward 1 /T^ as gain 
increases, while the two phugoid roots rendezvous and become a new 
second order; the short-period damping ratio increases somewhat, and 
the short-period undamped natural frequency remains essentially unchanged. 

The crossover region (frequency region near the intersection of the 
0 dB line and the Bode) compatible with good closed-loop response lies 
along the approximately —20 dB per decade slope starting at about l/Tg^ . 
Selecting a gain for this beam deviation loop involves a compromise between 
the dominant path mode (the phugoid) and the mode associated with the very 
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low frequency dipole pair. The higher the gain the less the effect of 
the dipole, yet the less well damped becomes the quadratic path mode. The 
gain selected is therefore made to compromise these factors. The gain 
value, K^, chosen results in a crossover frequency somewhat greater than 
0.2 rad/sec. The dominant mode with this gain is the path control quad- 
ratic (nee phugoid) with £p = 0.445 and = 0.465 rad/sec. The very- low- 
frequency root at 0.028 will be dominant in some degree of freedom, such 
as speed, although in the deviation response to a d command its effect 
will be partially removed by the lead at 1 /T^ = 0.035 . Thus we have 
achieved a system which exhibits a well damped, fairly rapid response 
path mode which is slightly contaminated by a very long time constant 
mode, together with a relatively high frequency, reasonably damped, 
short-period mode. 

Let us now consider imaginary modifications to these nominal plots. 

In the deviation loop closure, the range of permissible crossovers would 
be extended if the breakpoints at 1 /Tq^ and 1 /Tq 2 could be further separ- 
ated. Also the siggy segments of the Bode root- locus would then be moved 
down relative to the asymptotic plot, such that the attainable path mode 
damping ratio for a given gain, K^, would be increased. Unfortunately, 
with the feedbacks available in System C the maximum separation attainable 
is limited by l/T 01 and 1 /Tq 2 , and these can be approached in the attitude 
loop closure only at the expense of an underdamped short-period. However, 
this short-period deficiency which is developed in the inner-loop closure 
is partially made up in the outer-loop closure, where £ S p is slightly 

increased over £ * . If the airplane- alone had greater short-period 
S P 

damping, then a larger attitude loop gain could be used, thereby per- 
mitting an increased deviation loop bandwidth with the same damping ratio. 
Via this reasoning, a pitch rate damper would be useful even when the 
short-period damping ratio is large, as on this DC-8 example. 

When only the attitude and deviation are permitted as feedbacks, the 
gains given here are nearly optimum, in that the resulting responses are 
rapid and well damped and the change in these responses will be relatively 
insensitive to changes in many of the vehicle parameters. From the sensi- 
tivity viewpoint, the primary effects are those of 1 /T^ , which will modify 
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the very-low*- frequency closed- loop root; l/Tg^ "which profoundly affects 
the dominant quadratic path mode; and the short-period undamped natural 
frequency (<n S p) and damping ratio (£ ). The sensitivity of the closed- 

loop roots to other open- loop characteristics is not large. This is most 
easily appreciated by recognizing that the first-order sensitivities * 

(Ref. 2 )^ relating incremental changes in closed- loop roots to incremen- 
tal changes in open- loop gain or open- loop roots, are inversely propor- 
tional to the slopes on the Bode root-locus. As can be seen from Figs. .10 
and 1 1 , these slopes are quite large at the chosen gains . 

The closed- loop transfer functions for System C are given in Table 8. 
These data are used in Ref. 1 to compute the rms deviations due to u and 

O 

inputs. 

TABLE 8 

CLOSED-LOOP TRANSFER FUNCTIONS FOR SYSTEM C 

d £ +0. 566(0)(15.229)[0. 134 , 2.087] 

u g (0.028)(2.066)(15.228)[0.445 , o.465][o.206 , 2.039] 

d e +1.5(0. 110(15. 139)[0. 291 , 1.76] 

w g ( 0 . 028 ) ( 2 . 066 ) ( 1 5 • 228 ) [ 0 . 445 , 0.465] [0.206 , 2.039] 


d € 

^•command 


+2.0(0)(0.l3)(0.46)(l5.228)[0.l8 , 2.05] 
(0.028)(2.066)(15.228)[0A45 , 0.465] [0.206 , 2.039] 


u 0. 0373(0. 12)(1.35)(2.169)(15.228)[0. 176, 1 . 995 ] 

u g ( 0 . 028 )( 2 . 066 )( 15 . 228 )[ 0 . 445 , 0 . 465 ][ 0 . 206 , 2.039] 


u -0.136(2.051 )( 15 . 224 )[~ 0 . 17 , 0.442] [O.165 , 1 .96] 
w g (0.028) (2. 066) (15. 228) [0.445 , 0.465] [0.206 , 2.039] 

0.1 94(0) ( 4 . 03) (- 4 . 082) 

d command (0.028)(2.066)( 15 . 228 )[ 0.445 , 0 . 465 ] [0.206 , 2.039] 
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Typical transient responses for a deviation command, d c , and step 
Ug and Wg gust inputs are shown in Figs. 12 through l4. The step 
response to a d command is, as anticipated, dominated by the very- low- 
frequency dipole pair and the dominant path mode. Only a slight over- 
shoot is present with the gains selected and this response is, in 
general, entirely satisfactory. The short-period properties show up 
primarily in the elevator trace. The attitude primarily simulates the 
altitude rate, as can be seen by comparing the h and 0 traces. Finally, 
a small speed deviation occurs with a time constant given by the closed- 
loop mode, 1 /T^ , which approaches 1 /T^ . This is more graphically 
demonstrated on the speed responses to u g and w g step disturbances. 

System with Low Frequency Attitude Washout (System B) 

System B is very similar to System C. The difference lies in the 
attitude feedback at very low frequencies which is washed out rather than 
a pure gain. This has several effects. Statically, the need to fly with 
a slight deviation to offset any steady- state attitude is removed. Dynami- 
cally, the following of beam deviation commands and response to w gusts can 
be made somewhat better. 

The open- and closed-loop attitude dynamics are shown on Fig. 15 . When 
this is compared with Fig. 10, the high frequency characteristics are 
seen to be very similar, whereas the very-low-frequency properties are 
quite different. The open- loop dc gain is Eero, and the presence of 
the washout gives rise to a root at — l/T^o* For the same gain as used 
on the conventional system, the short-period frequency is essentially 
the same, the damping ratio is very slightly larger, and the phugoid 
roots are in closer proximity to 1 /Tq^ and 1 /Tq 2 . All of these features 
have favorable effects in the outer loop. 

The path deviation loop characteristics, with the attitude inner 
loop closed, are given in Fig. 16. By comparison with the conventional 
system, the low-frequency properties are close to those of a K/s system 
over a very wide frequency band. Consequently, the system gain can be 
set at almost any value in this range and result in a well- shaped response 
to a d c command. The time scale of the response is, of course, scaled by 
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Figure 16. System Survey for Beam Deviation Loop Closure for System B 


the crossover frequency. The same deviation- loop gain is used for this 
system as for the conventional. The closed- loop roots resulting are 
quite similar to those of System C, although the undamped natural fre- 
quency, o^, and damping ratio, £p* of the dominant mode are slightly less. 
This is to be expected as a consequence of the introduction of the 
washout. The closed- loop transfer functions are given in Table 9. 


TABLE 9 

CLOSED- LOOP TRANSFER FUNCTIONS FOR SYSTEM B 

+0.566(0)(0.019)( 15.23) [0. 1^7 , 2.099] 

U g (0.059)(0.07)(2.065)(15.229)[0.421)- , 0.4l5][0.2l8 , 2.06] 

a +i .5(0.015 ) (0.095 )( 15. i 4 ) [0.509 , 1.786] 

€ _ 

w g ( 0. 039) (0.07)(2. 065) (15. 229) [0.4214. , 0.415] [0.218 , 2.06] 



^■command 


+2. 0(0) (0.018) (0.109) ( 0 . 489 ) (15. 229) [0.20 , 2.07] 

( 0.039) (0.07) (2. 065 )( 15.229) [0.424 , 0.415] [0.218 , 2.06] 


0 . 0373 ( 1 . 366 )( 2. 1 67) ( 15. 229) [0.609 , O.O97] [O.191 , 2 . 013 ] 
u g ( 0.039)(0. 07) ( 2.065 )( 15.229) [0.424 , 0.4i5][0.2i8 , 2.06] 

-0. 136(0. 078)(2.05)(i5.225)[-0.24l , 0.444] [0. 180 , 1.978] 
w g (0. 039)(o. 07)(2. 065 )( 15. 229) [0.424 , 0.415] [0.218 , 2.06] 

0.1 94(0) (0.08) ( 4 . 03) (- 4 . 082) 

^command ( 0 . 039 )( 0 . 07 )( 2 . 065 )(l 5 - 229 )[ 0.424 , 0 . 4 l 5 ][ 0 . 2 i 8 , 2.06] 


The system transient response characteristics are generally similar 
to those of System C, with slight variations as expected by the qualita- 
tive system differences. The responses to a step d c command are almost 
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exactly identical. Some differences would occur, however, if the input 
were a higher order function, such as a ramp; then the response of 
System B would "be superior. For the w step inputs. System B is markedly 

o 

"better than the conventional System C in terms of beam deviation. The 
other degrees of freedom are generally similar, as would again be expected. 
Finally, with a u g step disturbance. System B is again superior in the d 
response, with the other degrees of freedom again very similar. Thus, in 
every respect but path deviation loop bandwidth, the modified system (B) 
is better than the conventional controller (C). 

Advanced System (System A) 

The advanced system is analytically more complex than the others 
considered because of the more complicated equalization and because the 
use of a K^h signal for path damping requires an additional loop closure. 
For simplicity of explanation, however, we shall use a K^d signal rather 
than a K^h signal. With this simplifying assumption we can again deal 
with an attitude and path deviation set of loop closures for the Bode 
plots. However, it is noted that for the closed- loop analog computer 
time responses, the actual K^h signal was used, rather than the simplified 
K£d signals. 

In the attitude loop, the attitude washout is made somewhat less than 
the short-period undamped natural frequency so as to assure nearly pure 
gain attitude feedback at short-period frequency. The addition of the 
attitude rate signal to the washed out attitude creates a net equalization 
on 0 given by 


where 


K 0 s 


s + 1 /T wo 


+ K a s 


k 0 s(s + 1/afc) 
S + 1/Two 


j_ _ fe , 

t E K 0 t wo 


(41) 
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The location of the lead equalization breakpoint, 1 /T e , is the primary 
means to adjust the short-period damping ratio. If this breakpoint is 
placed somewhat greater than <n , a long stretch of —20 d B per decade 
slope will be established between 1 /Tg and the actuator breakpoint at 
l/T a . Gain crossover anywhere in this stretch will result in reasonable 
damping of the short-period mode, while still permitting reasonably large 
amplitude ratio values at mid- frequencies . This is illustrated in the 
attitude control jo> Bode and Bode root-locus in Fig. 17. There, it is 
seen that with an increase in gain the phugoid and short-period damping 
ratios, and the short-period undamped natural frequency are all increased, 
while the phugoid undamped natural frequency is decreased. The physical 
explanation for the decrease of phugoid frequency with System A is the sig- 
nificantly lower amplitude ratio of the e feedback at phugoid frequency 
(compared to System B), which is due to the difference in the washout time 
constants . 

Consider now the beam deviation closure. Here, the equalization, has 
both a rate and integral term in addition to the glide- slope-beam noise filte 

K 3 + K d s 

+ K;s 

s(T f s + 1 ) 

K d [ s 3 + ( 1 /T f ) s 2 + ( Kfl/Tf K£ ) s + K 5 /T f K d ] 
s(s + l/T f ) 

K d ( s + 1 /T dl )(s + 1 /T d2 )( s + 1 /T d3 ) 
s(s + 1 /T f ) 

fd J_ ^ % _1_ ^ _1_ 

K a ' % = ’ T *5 = T f 

2 

At very low frequencies, the open-loop system looks like K/s because of 
the integral term in the deviation controller. The rest of the equaliza- 
tion is needed to adjust the amplitude ratio to approximate a —20 dB per 
decade slope in the desired region of crossover. While has three time 
constants, 1 /T d ^ is inherently close to the beam noise filter breakpoint. 


i.e.. 


where 


1 

Td! 
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so the best that can be achieved is to create a dipole pair near 1 /T f 
giving a small phase lead in the region of the filter. The principal 
adjustments are then 1 /T^ and l/T^. Wh e & these breakpoints are posi- 
tioned relative to the breakpoint, as shown in Fig. 18, the low- 
frequency amplitude ratio attainable with a given zero dB line is made 
as large as possible. However, note that the phase in the low-to-medium 
frequency region must not be permitted to exceed — 180°. This is impor- 
tant for avoiding a low-frequency, low-amplitude oscillation due to any 
system threshold or hysteresis characteristics (either of these will 
result in reduced gain at very low amplitudes which could set the condi- 
tions for a limit cycle at a frequency where the phase angle reaches —180°). 
If the loops are closed with the gains noted, then a crossover frequency 
corresponding to a maximum attainable phase margin of about 32° is achieved. 
The resulting closed- loop characteristics are greatly superior to those of 
either the conventional or modified system in bandwidth, path following, 
dominant and short-period mode characteristics, etc. The closed- loop 
transfer functions are given in Table 10. 

The general superiority of System A is also exhibited in the transient 
response comparisons of Figs. 12 through '\k. Deviation responses for both 
command and disturbance inputs are all superior, as are all the degrees of 
freedom which most strongly reflect short-period properties. A summary 
table of the various control system constants is given in Table 11 for 
easy reference. 

The speed responses are very similar for all three systems. This 
reflects the lack of a speed control loop. The basic time-constant of 
the speed deviation is limited by l/T^. In terms of approximate factors, 
this is given by 


1 


-x u + (Xa-g) 


a* 0 -Z5 e M u /M 5e Z u ) 

^ 0-Z 5e Hx/M &e Za) 


(* 2 ) 


When a speed-control system is considered as a stability augmenter, an a 
or u to 5-5 feedback will modify X u or X^ and thus can be used to change 
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1 / Th -j • Then the deviation loop closure will result in a larger value of 
1 /T]^ , with a concomitant improvement in the closed-loop system speed 
response in all three systems. 

TABLE 10 

CLOSED -LOOP TRANSFER FUNCTIONS FOR SYSTEM A 

d e _ +o. 566(0) (o)(o. 17*0 (12.918) [0.767 , 2.215] 

u g (0.036)(0.i23)(0.582)(2.462)(i3.232)[0.657, O.699H0.673 , 1.428] 

_ +1 -3(Q)( 13.75) [0.464 , 0.105] [0.936 , 2.018] 

w g (0.036)(0.123)(0.582)(2.462)(13.232)[0.657 , 0.699] [0.673 , 1.428] 

d e _ +2.0(0) 2 (0.044)(0.733)(2.0)(13.22)[0.43 , 1.46] 

^command (0.036) (0. 1 23)(0.582)(2.462)( 1 3.232) [0. 657 , 0.699][0.673 , 1 .428] 

-0. 0373(0. i36)(i.596)(2.777)(i3. 261 )[0.5 , 0.276H0.58, 1.918] 
u g ( 0.036) ( 0.1 23) (0. 582 )( 2.462) ( 13.232) [0.657 , 0.699] [0.673 , 1 .428] 

-0. 136(0. 153)(0.215)(2.409)(i 3.262)[-0. 082, 1 .023] [0.872 , 1.492] 
w g ( 0.036) ( 0.1 23) (O.582) (2. 462) ( 13. 232) [0.657 , 0.699] [0.673 > 1.428] 

0.3272(0) (O.O89) (0.7) (-4. 082) (4. 03) 

^command (0.036)(0.123)(0.582)(2.462)( 1 3.232) [0.657 , 0. 699] [0.673 , 1.428] 
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TABLE 11 


SUMMARY OF CONTROL SYSTEM CONSTANTS 



A 

B 

C 

K e 

-2. 

-3.652 

-3.652 

Kq (sec) 

—2. 

0 

0 

(rad-sec/ft) 

-0.0256 

0 

0 

(rad/ft) 

-0.00867 

-0.00514 

-0.00514 

Kj (rad/ft-sec) 

-0.000768 

0 

0 

l/T wo ( 1 /sec) 

0.7 

0.08 

0 

1 /T a ( 1 /sec) 

15. 

15. 

15. 

1 /T f ( 1 /sec) 

2. 

2. 

2. 
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